clearvars
clc

% Produce Figure 8

load('ResultsDeltaTable')
deltatem = 0.990721; % This is the delta from the baseline calibration
Dinf     = sum(USInfhu(:,2:11),2)'-sum(USInfhu(27,2:11),2)'+(deltatem-delvec)*1000;
Dune     = sum(CUSUnehu(:,2:11),2)'-sum(CUSUnehu(27,2:11),2)';

figure(8)
plot(Dune,Dinf)
title('Unemployment vs. inflation ``tradeoff"','Interpreter','Latex')
xlabel('U.S. 2001-2010 unemployment relative to baseline','Interpreter','Latex')
ylabel('U.S. 2001-2010 inflation relative to baseline','Interpreter','Latex')
axis([-6 2 -5 50])
